Stability of Bose Einstein condensates of hot magnons in YIG 
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We investigate the stability of the recently discovered room temperature Bose-Einstein condensate 
(BEC) of magnons in Ytrrium Iron Garnet (YIG) films. We show that magnon-magnon interactions 
depend strongly on the external field orientation, and that the BEC in current experiments is actually 
metastable - it only survives because of finite size effects, and because the BEC density is very low. 
On the other hand a strong field applied perpendicular to the sample plane leads to a repulsive 
magnon-magnon interaction; we predict that a high-density magnon BEC can then be formed in 
this perpendicular field geometry. 

PACS numbers: 



In a remarkable very recent discovery [l[ , a Room Tem- 
perature Bose-Einstein condensate (BEC) of magnon ex- 
citations was stabilized for a period of roughly 1 fis in 
a thin slab of the well-known insulating magnet Yttrium 
Iron Garnet (YIG). The density of magnons was quite low 
(n ~ 10~ 4 per lattice site), and the density n of the BEC 
was apparently unknown, but n /n <C 1. This result can 
be understood naively in terms of a weakly-interacting 
dilute gas of bosons, provided that (i) one assumes that 
the number of magnons is conserved, so their chemical 
potential may be non-zero, and (ii) the interactions be- 
tween them are repulsive (attractive interactions favour 
depiction of the BEC, causing a negative compressibility 
and instability of the BEC [2|]). In the experiment, the 
magnon dispersion was controlled both by the sample ge- 
ometry and external magnetic field, in such a way that 
magnon-magnon collisions conserved magnon number - 
the decay of the BEC was then attributed to spin-phonon 
couplings [l[ . The magnons in the experiment had rather 
long wavelengths, of order a /jm - hitherto such excita- 
tions have been treated entirely classically 

This experiment raises a number of important ques- 
tions, not least of which concern the kind of supcrfluid 
properties possessed by a BEC of such unusual objects. 
However a much more basic question about the stability 
of the system must first be answered. In fact we find the 
rather startling result that under the conditions of the 
experiments reported so far, these interactions were ac- 
tually attractive, ie., the BEC ought to be unstable! We 
shall see that because of the particular geometry used, 
the BEC is actually metastable to thermally activated 
or tunneling decay, and that it only survives because its 
density is low - above a critical density, given below, it is 
absolutely unstable. However we also find that by chang- 
ing the field configuration in the system one can make the 
interactions repulsive, and the BEC should then stabilize 
at a much higher density - opening the way to much more 
interesting experiments. 

YIG is one of the best characterized of all insulat- 



ing magnets It is cubic, with lattice constant a = 
12.376A, ordering ferrimagnctically below T c = 560A'. 
At room temperature the long-wavelength properties can 
be understood using a Hamiltonian with ferromagnetic 
exchange interactions between effective 'block spins' Sj, 



one per unit cell, whose magnitude Sj 



a 3 M a />y, 



with 7 = g e fiB, is defined by the experimental saturated 
magnetisation density M s (M s m 140 G at room temper- 
ature; with g e w 2, one has Sj ~ 14.3), along with dipolar 
couplings between these; the resulting lattice Hamilto- 
nian takes the form: 



H 



(i) 



where the sums i, j are taken over lattice sites at po- 
sitions R;, etc., 5 denotes nearest- neighbor spins, = 
(R,— Rj)/a , and = r^-/|rjj|. The nearest-neighbour 
dipolar interaction U d = r ) 2 /a Q k 1.3 x 10~ 3 A'. The 
isotropic exchange J a is determined experimentally from 
J = 2SJ a f=a 0.83 x 10~ 28 erg cm 2 at room tem- 
perature. One then has J « 1.37 K, and Ud/J ~ 
0.95 x 10~ 3 . 

In what follows we set up a theoretical description of 
the BEC, taking into account the external field, dipolar 
and exchange interactions, and boundary conditions in 
the finite geometry. We evaluate the interactions and 
the BEC stability for 2 different field configurations; the 
general picture then becomes clear. 

(i) In-plane field: All experiments so far have had 
H D in the slab plane. The combination of exchange, 
dipolar, and Zceman couplings then leads to a magnon 
spectrum u q shown in Fig[TJ in which the competition 
between dipolar and exchange interactions leads to a 
finite-g minimum in o; q at a wave- vector Q ~ 1 /d, where 
d is the slab thickness. To completely specify w q and 
the inter-magnon interactions one needs boundary con- 



2 



ditions, which can involve partial pinning of the surface 
spins [a, & Q ■ Demokritov et al [l[ assume free surface 
spins, implying that (i) \(dM(r)/dr) ■ n s \ = when r is 
at the surface; here n s is the normal to the surface, and 
(ii) that the allowed momenta along z (see Fig[TJ inset) 
are q± = n±n/d, leading to different magnon branches 
labelled by n±. For now we assume a continuous in-plane 
momentum, and later discuss the effect of in-plane quan- 
tization; and we assume n± = 0, taking the lowest-energy 
magnon branch. 
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FIG. 1: (color online). The magnon spectrum, EqO for 
d = 2, 5 and 10 /im at H = 700 Gauss. The inset shows the 
sample geometry. 

The above assumptions yield an oj q which agrees with 
experiment [H, and which for H G || x takes the form Q 



hio q = [( 7 ff t + Jq 2 )("/Hi + Jq 2 + tiUMFq)] 



1/2 



(2) 



here Hi = H — 4irN x M s is the internal field, N x the 
demagnetization factor (for a slab in the xy-plane, N x — 
N y = 0, N z = 1), tvjj M = 4ir~/M s (for YIG, huo M ~ 
0.236 K at room temperature), and hoj H = ■yHi. The 
form of the dimensionless function T q depends on the 
direction of q. In the important case when q is parallel 
to H a , ie., along x, one has 



T 

J n 



'1-e 



/(q x d), 



(3) 



In the experiment [l| magnons are argued to condense 
at the minima q x = Q; when d — 5 [im, and H a = 
700 Gauss, one has \Q\ ps 5.5 x 10 4 cmT 1 . 

We now set up a theoretical description of the BEC, 
including all 4-magnon scattering processes (3-magnon 
scattering is excluded by the kinematics), using a gener- 
alized Bogoliubov quasi-average technique 0] to incor- 
porate the BEC. Defining magnon operators 6 q ,6q, a 
magnon BEC at q = Q, with N a condensed magnons, 
has quasi-averages 



< b ±Q >=< bl Q >= s/N^/2, 



(4) 



corresponding to a condensate wave- function ^q(u) oc 
cos(Qy). More generally v E , q(y) is multiplied by a phase 
factor e I< ^ r, '- ) / ? ', which is crucial to the BEC dynamics, 
but not necessary for a stability analysis of the BEC. 

We make a Holstein-Primakoff magnon expansion Q 
up to 4th-order in magnon operators, including contribu- 
tions from both (2 in - 2 out) and (3 in - 1 out) magnon 
scattering processes (we ignore multiple-scattering 
contributions here, which are ~ 0(Ud/J o ) ~ 10 -3 rela- 
tive to the leading terms). Then, taking quasi-averages, 
we can write the Hamiltonian in the form 



VI 



(5) 



where the interaction term V[ nt takes the form 



n (T 
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^ [bq+pbn-p + fc-Q-p^-Q+p + 
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Here To and T$ are the four-magnon scattering ampli- 
tudes between states (Q, Q) — » (Q, Q) and (Q, — Q) — > 
(Q, — Q) respectively. For the sample geometry in Fig[TJ 
with H Q || x, these scattering amplitudes arc found to be 



To 
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[(ai - u z )T Q - 2a 2 (l - Fiq) 
JQ 2 
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s ■ [ai - 4a 2 ] ; (7) 
[(ai - a 2 )(l - T 2 q) - {pt\ - a 3 )T Q ] 



JQ 2 , 
+ -q-[<xi 



2a 2 



(8) 



with a± = Uq + AuqVq + Vq, a 2 = 2uqVq and a 3 = 

1 /2 

3uqVq(uq + Vq), where {u q ,v q } = [(A q ± frw q )/2fruj q \ 
and ^4 g and B g are given by A q = ■yHi + Jq 2 + 2it~fM s F q 
and B q = —ir^AI s J- q respectively. Higher-order multiple- 
scattering contributions to ro,rs are ~ 0(Ud/J ) rela- 
tive to the leading terms given here. 

We plot the combinations 2Fo±rs in Figures [2] and [5] 
Both these amplitudes are sensitive to the external field 
and the film thickness. The first amplitude, 2To + Tg, 
becomes negative in the entire region of fields at d < 
d c « 2(J/huj M ) l/2 « 0.032 fim. The second amplitude is 
positive at d < d c . 

Near the energy minimum (when \p\ << |Q|), one has 
WQ+p ~ LOQ-p, and the Bogoliuibov transformation is 



3 



straightforward because 7i is symmetric when p <-> — p 
and Q <-> — Q. The spectrum thus has 4 branches, with 
excited quasiparticle energies 



tpo = \Z^q(p)[^q(p) +?^o(2r +77r s )] 
where 77 = ±1, and £Iq(p) = h(aj Q+p —uiq). 

2r +r s 
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FIG. 2: (color online). The effective amplitude 2Fo + Ys (in 
Kelvins) as a function of field H at different values of the 
film thickness d. Dashed lines: d = 2 and f /im. Solid lines: 
d = O.Of , 0.05 and 5 fj,m. 

The inset in Fig[3] shows the "phase diagram" of the 
quasi-two dimensional YIG for different values of d and 
Hg . One immediately sees a paradox: the amplitudes are 
never both positive, so bulk BEC should not exist - yet 
BEC has been observed [lj in samples with an in-plane 
field H% ~ 0.07 T, in which d ~ 5 pm. 

The paradox is resolved by noting that in a finite ge- 
ometry, the energy gap from the condensate to excited 
modes can be larger than the scattering amplitude, lead- 
ing to a potential barrier to decay. For weakly-interacting 
Bose gases this yields 11| an upper critical number n^ r of 
condensate particles, above which the barrier disappears; 
one has ot n c J jl = k, where a is the s-wave scattering 
length, and l a the characteristic size of the BEC wave 
function. The constant k ~ 0(1), and depends on the 
sample geometry. 

In the present case we can write the critical density 



as \2To~TsK 



where 



is the minimum 



quasiparticle energy in the presence of the BEC; below 
this critical density the BEC is metastable to tunneling or 
thermal activation. If the BEC were to spread through 
the entire slab, then n^ r \2To — Ts\ = Qq(tt/L), where 
the length L depends on the direction of the field rela- 
tive to the slab axes. Taking this result literally for the 
experiment [l[ , with a slab measuring 20x2 mm 2 in the 
plane, one finds 10~ 8 < n c J < 10~ 6 (for fields along the 



long and short sides of the slab respectively). However 
this result is certainly too low, since it assumes a per- 
fectly uniform BEC - in reality disorder and edge effects 
will smear the magnon spectrum and restrict the size 
of the BEC, and a more realistic estimate for L is then 
L eff = (LxLyLz) 1 / 3 ~ 0.6 mm. This yields nf ~ 10~ 5 
for the experiment. Thus we conclude that for in-plane 
fields, it will be impossible to raise n above this value; 
this could be checked experimentally (eg., by increasing 
the pumping rate). 
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FIG. 3: (color online). The effective amplitude 2Fn — Ys (in 
Kelvins) as a function of external field H at different values 
of the film thickness d. Dashed lines: d = 2 and 10 fim. Solid 
line: d — 0.01, 0.05 and 5 \xm. The inset shows the "phase 
diagram" of YIG in the (H , d) plane. The amplitude 2ro+rs 
is positive everywhere to the right of the solid line (the region 
©). The amplitude 2Fn — Ys is negative everywhere to the 
right of the dashed line (the region 0). 

(ii) Perpendicular field: Now the results are very 
different. Consider first an infinite thin slab, which is 
simple to analyze. The competition between the external 
field H = zH a and the demagnetization field (which 
favours in-plane magnetisation) gradually pulls the spins 
out of the plane; below a critical field H c = hwu/l ~ 
1760 G, the Free Energy is degenerate with respect to 
rotation around z and so the magnons are gapless, but 
at H c they align with H Q and a gap fku = ftw q= Q 
" II, — huiM opens up. The minimum in the magnon 
spectrum is always at q = (FigHJ). 

The inter-magnon scattering amplitude T is now al- 
ways positive; neglecting a very small exchange contri- 
bution one finds 



r(q) = (^ M /45)[l-(l-^)/2] 

-» (tiu M /4S)[l-q n d/4: + 0{q 2 )) 



(10) 



where q\\ is the momentum in the xy plane, and T q takes 
the form §3§ but with q x — > q\\ . This radically changes 
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the situation - now a BEC is stable with no restriction 
on the condensate density. There are however restric- 
tions on H ; when 2200 G < H* < 3500 G the system 
has a "kinetic instability" [T2|, in which the pumping of 
the magnons at one frequency destabilizes the magnon 
distribution, along with strong microwave emission. 
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FIG. 4: (color online). The magnon spectrum for an infinite 
slab of YIG, assuming free surface spin boundary conditions, 
with magnetisation polarised perpendicular to the plane (see 
inset). The spectra are shown for HI — H c ~ 1760 G where 
the spectrum is gapless, and for H* = 2000 G. 

In a real finite sample things are much more compli- 
cated. Even without surface anisotropy the spins near the 
slab edge are put out of alignment with the bulk spins by 
edge demagnetisation fields; and surface anisotropy does 
the same to spins on the slab faces. However in the cen- 
tral region of the sample, at distances further from the 
surface than the exchange length, the spectrum returns 
to the infinite plane form. For fields well above H c , eg., 
for H* ~ 2000 G, all of the spins will be aligned along z, 
and pop will then be valid everywhere. 

In this case we have the striking result that a BEC 
of pumped magnons should be possible with densities 
much higher than present. To give an upper bound is 
complicated since the problem then becomes essentially 
non-perturbativc (similar to, eg., liquid He), beyond 
the range of higher-order magnon expansions. However 
there appears to be no obstacle in principle to raising 
n /n ~ 0(1). At present the highest achievable density 
is probably limited by experimental pumping strengths 
rather than any fundamental restrictions. Such a high- 
density BEC existing at room temperature would be ex- 
tremely interesting, and certainly possess unusual mag- 
netic properties. 

Remarks: The 2 cases studied above are actually lim- 
iting cases of a more general situation in which one can 
manipulate the inter-magnon interactions by varying the 
field direction and strength, and vary the upper criti- 
cal density for BEC formation by changing the sample 



geometry. Thus the analysis can be easily generalized 
to long 'magnetic wires' or whiskers, and we also ex- 
pect that BEC will be stabilized there when the external 
field is perpendicular to the sample axis, but unstable or 
metastable when the field is parallel to the sample axis. 
Further details of the various possible cases will be pub- 
lished elsewhere. 
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